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Abstract. After formulating the frequency-domain Maxwell equations for a homogeneous, 
linear, bianisotropic material occupying a bounded region, we found that the axionic piece 
vanishes from both the differential equations valid in the region and the boundary conditions, 
thereby vindicating the Post constraint. Our analysis indicates that characteristic effects that 
may be observed experimentally with magnetoelectric materials are not the consequences of 
the axionic piece but of an admittance that describes surface states. 
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1 Introduction 

Do axions [1] exist? We, the authors of this communication, lack the capabilities to answer 
that question. But we do note several negative answers to that question have emerged 
from experiments [2-5]. Does the “axionic piece” [6] of the linear constitutive dyadic exist 
in modern classical electromagnetic theory? Again, we are not in a position to answer 
this question. But one of us summarized the developments before 2004 to state that the 
recognizable existence of the axionic piece is ruled out by modern classical electromagnetic 
theory [7]. The filtering out of the axionic piece (denoted by T in this communication) by 
the Maxwell equations is enshrined as the Post constraint T = 0 [8,9], having been delivered 
by Post [10] more than half a century ago. 

The validity of the Post constraint has been vigorously debated. A simple boundary- 
value problem has been formulated to show that T is measurable [6]. Multipole symmetry 
arguments have been derived to rule against as well as for the Post constraint [11-13]. Data 
from scattering experiments conducted on a single “Tellegen particle” were produced to 
show that T exists [14], but those data—while sufficient to establish the magnetoelectric 
phenomenon (which is not in doubt)—are insufficient to invalidate the Post constraint. Al¬ 
though suspensions of Tellegen particles have also been shown to exhibit the magnetoelectric 
phenomenon [15,16], no data has been presented to invalidate the Post constraint. 

There is one item of experimental evidence against the Post constraint, as recounted in 
great detail by Hehl et al. [17] in 2007. This evidence emerging from (i) certain assumptions of 
symmetry of the linear constitutive dyadic that lead to the description of the magnetoelectric 
properties of Cr 203 in terms of just two scalars; (ii) DC measurements reported in 1994 to 
obtain the sign and the magnitude of one of the two magnetoelectric scalars as well as 
the magnitude of other magnetoelectric scalar of Cr 203 on a disk-shaped sample and a 
rectangular-solid sample [18], respectively; (hi) the sign of the second magnetoelectric scalar 
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of Cr 203 measured at 10 kHz on a spherical sample reported in 1961 [19]; and (iv) the 
experimentally supported approximation that the quasistatic permeability of Cr 203 is the 
same as that of vacuum [17,20], enabled them to deduce non-zero values of H/. Although 
independent conhrmation of this experimental evidence does not exist to our knowledge, the 
veracity of measured data from careful experiments should not be doubted. 

But the interpretation of measured data in order to obtain derivative quantities of greater 
signihcance is not necessarily unimpeachable. It is our contention here that data from 
electromagnetic experiments on samples of magnetoelectric materials (such as Cr 203 ) do 
not yield evidence in favor of non-zero H/. Instead, the effects attributed to a non-zero tl', 
a macroscopic quantity that supposedly exists throughout a sample of a magnetoelectric 
material, arise from surface states that exist due to the abrupt cessation of microscopic 
morphology at the boundary of that sample. 

The remainder of this communication is planned as follows: Section 2 is devoted to a 
mathematical exposition of the Post constraint being a natural outcome of the application 
of Maxwell equations to linear materials. As T does not appear in the Maxwell equations, it 
is necessarily absent from the boundary conditions derived from those differential equations 
[21]. Section 3 presents the macroscopic characterization of surface states in terms of a 
surface charge density and a surface current density mediated by an admittance. Section 4 
contains the solution of a boundary-value problem to show that certain characteristic effects 
attributed to the axionic piece are the consequences of that admittance. This communication 
concludes with some remarks in Sec. 5. Vectors are underlined, 3x3 dyadics [22] are double 
underlined, 0 is the null vector, / is the 3x3 identity dyadic, and 0 is the 3x3 null dyadic, 
£o is the permittivity of vacuum, and /Iq is the permeability of vacuum. 


2 Natural emergence of the Post constraint 

Suppose that all space is divided into two mutually disjoint regions Vout and as well as 
the surface S that separates Vout and Vin- We formally distinguish between Vouti Wn, and S. 
A differential equation equation holding in Vout (or Vin) would have to be solved to satisfy 
boundary values of dependent variables prescribed on S. 

Suppose that the region Vout is vacuous, and the region Vin is occupied by a homoge¬ 
neous, linear, bianisotropic material. With an exp{—iut) dependence on time t with angular 
frequency u and i = \/— 1 , the frequency-domain constitutive relations of free space are 
specihed as 


D{r, u) = eoE{r, u ), H{r, u) = VoBiji, w), re Vout , (1) 

with Uq = The frequency-domain constitutive relations of the bianisotropic material 

are specihed as 

= £(i^) • E_{Li + a(i^) • cn) -|- T(ci;)H(r, oj) 

= /^(<^) ' E{r,u) + ^(cu) • B_{r,u) — d^(cu)E(r, cu) 


r G Vi, 


( 2 ) 
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where 


Trace 


«(ci;) — I3{uj) 


= 0 


( 3 ) 


The 3x3 permittivity dyadic e, the 3x3 impermeability dyadic the two 3x3 magneto¬ 
electric dyadics a and /?, and the axionic piece T can be combined into a single 6x6 linear 
constitutive dyadic [22], but we do not use that compressed notation here. When writing 
the constitutive relations (2) for D and ^ in terms of ^ and 5, we have kept in mind that 
^ and 5 are primitive fields because of their underlying microscopic existence while D_ and 
H are induction fields that arise at the macroscopic level but do not exist at the micro¬ 
scopic level. The non-existence of induction helds at the microscopic level is a cornerstone of 
modern physics—in contrast to a heterodox formulation of Hehl and Obukhov [23] wherein 
the merely convenient but inessential D and ^ at the macroscopic level have microscopic 
counterparts, thereby obscuring the distinction between the microscopic and macroscopic 
levels. 

Let us now apply the frequency-domain macroscopic Maxwell equations 


V • B(r, cu) = 0 

V X E(r, oj) — = 0 

V ’ D_{r,u) = p{r,u) 

V X H_{r, u) -I- iu}D{r, uj) = J(r, uj) 

to both Vout and Vin- After substituting Eqs. (1) in Eqs. (4) we get 

V • 5(r, (u) = 0 

V X E{r, u) — cu) = 0 

> , r G 

ffoV ' E{r,uj) = p{r,u) 

1^0V X ^{r,uj) + iuEoKiE^) = ULi^) 


( 4 ) 


( 5 ) 


Likewise, after substituting Eqs. (2) in Eqs. (4) we get 


V • S(r, cu) = 0 

V X E(r,a;) — cu) = 0 

V • \g{uj) • E(r, uj) -f a{uj) • 5(r, cu) -|- T(a;)5(r, cu)] = p(r, uj) 


> , r G Vjn , 


V X 


(3{uj) ’ E{r, uj) + u{uj) ’ B_{r,uj) —'^{uj)E{r,uj) 

+iu [£(0;) • E_{r, u) + a{uj) • 5(r, cu) -|- T(ci;)5(r, a;)] = J(r, uj) 


( 6 ) 
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which can be first rearranged as 

V • R{r,uj) = 0 

V X E(r, u) — iujB_{r, = 0 

V • [£( 0 ;) • + fl/(ci;)V • R{r,uj) = p{r,uj) 


> , r e Vi, 


V X 


/3{uj) • :S(r, cu) + ^(cu) 'B_{r,uj) — \I/(a;) [V x E(r, cu) — ia;5(r, cu)] 
+iu [ 1 .( 0 ;) • ^(r, u) + a(a;) • 5(r, a;)] = J(r, a;) 


and then simplified to 

V • B(r, w) = 0 

V X E(r, a;) — iuj^{r, 1 ^) = Q 

V • [e(ci;) • E{r, uj) + a{oj) • 5(r, w)] = p{r, u) 

V X 


( 7 ) 


r e Vi, 


( 8 ) 


13(oo) • E(r, 00 ) + pipj) • 5(r, oj) 

+iu [e(a;) ' E{r,u) + a{u}) ' B_{r,u})] = J_{r,oj) 

Most importantly, does not appear in the Maxwell equations applied to Vin when the 
merely convenient induction fields have been translated into primitive fields. 

When solving an electromagnetic boundary-value problem, it is common to use the fol¬ 
lowing boundary conditions derived from Eqs. (4): 

fkir,) • = 0 

ikiLs) X [E^tiLs,^) - Ei^iLs,^)] = Q 

niLs) • = PsiLs,^) 

n(rj X 


hs e 5, 


(9) 


with the unit vector h(r^) at r^ G iS pointing into Vout- The quantities Ps and are the 
surface charge density and the surface current density, respectively. The subscripts in and 
out indicate that the fields are being evaluated on the interior and the exterior sides of S. 
The boundary conditions (9) are derived by integrating Eqs. (4) over pillboxes and closed 
contours straddling S, as appropriate. Substituting the constitutive relations (1) and (2) 
in Eqs. (9) and enforcing the charge neutrality and current neutrality of S, we get as the 
boundary conditions 

ZLiLs) • [Bout{rs,uj) - = 0 

n{r^) X [E^tiLs,(^) - E^^{r^,u)] = 0 

^oE^uALs, • E^^{r^, u) 

- T(a;)5^„(r„a;) = 0 
VoBouAl^, uj) - ^{uj) • a;) 

-k{uj) -Bi^ir^u) + T(u;)E.„(r,,u;) = 0 


mr. 


n(rj X 


, 


hs e 5, 


( 10 ) 
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which clearly contain \1/. Thus, one could conclude that although T has vanished from 
the Maxwell equations applied to Vin, it still makes its presence felt through the boundary 
conditions (10). 

That conclusion would be erroneous [21], As stated elsewhere [24], “the actual [emphasis 
in the original] issue is the correct physical formulation of boundary-value problems.” The 
boundary conditions (9) are not sacrosanct in electromagnetics, but their emergence from 
the Maxwell equations is. Using the same pillboxes and closed contours as in the derivation 
procedure for the previous paragraph, but using 


V • B(r, cu) = 0 

V X E{r, oj) — cu) = 0 

V • |l(z:, w) ’ ’ B_{r_,uj)\ = p{r,u) 


V X 


(3{r,u)) • E{r,u) +p,(r,uj) • 

+iu [l(r, a;) • :U(r, cu)-|-«(r, a;) ' = J_{r,uj) 


with 




f \ 

f 0 

r\ a{r,u}) = 


[ i(^) ^ 


Jo ^ J 

uoL 

1 = ■! 



r G 


Uout 


instead of Eqs. (4) because T vanishes from Eqs. (6) to deliver Eqs. (8), we get 


UlLs) • [RoutiZs,^) - = 0 

n{r^) X [E_^tiZs,uj) - E^Rr^,uj)] =0 

UZs) • uj) - g{uj) • E^{r^, uj) - • RiRZs, ^^)] = 0 ^ 

UZs) X l^oRoutiZs, • RiniZs, • RiniZs, Uj) = 0 


( 11 ) 


( 12 ) 


(13) 


as the boundary conditions. These boundary conditions do not contain T, in contrast to 
the argument raised elsewhere [25] based on the applicability of the Maxwell equations on 
S with T multiplied by the unit step function with zero value for r G Vout- 

What is the reason for the differences between the boundary conditions (10) and (13)? 
Very simply, the boundary conditions (10) were derived from Eqs. (4) without considering 
that the redundancies possibly contained in constitutive equations could be filtered out by 
the Maxwell equations when applied separately to Vm and Vout- When the redundancy in 
the form of 4/ was filtered out by the Maxwell equations on application to Vm, we obtained 
the boundary conditions (13) [21]. Thus, 4/ does not have a recognizable existence in either 
Vin or Vout and it does not enter the boundary conditions. Ergo, the Post constraint 


T(a;) = 0 


( 14 ) 


emerges as a natural outcome of the application of Maxwell equations to linear materials. 
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Let us stress here that the standard boundary conditions in macroscopic electromagnetics 
were not changed—in order to exclude d'. That constitutive parameter was absent in the 
applicable differential equations (i.e., the Maxwell equations containing the primitive but 
not the induction fields) and therefore did not appear in the boundary conditions derived 
from those differential equations [21], Any suggestion that the foregoing procedure “actually 
changes [emphasis in the original] the physical laws, namely, the Maxwell equations at the 
interface between the two [homogeneous] media” [26] appears strange to us because any 
differential equation is held to be valid in a region but not on the boundary of that region. 

Remark 1: If Vin were to be occupied by a homogeneous, linear, biisotropic material, 
then e = eJ, z/ = i/J, a = al_ = P, and T = 0. In other words, a biisotropic material must 
be Lorentz reciprocal and nonreciprocal biisotropy does not have a recognizable existence in 
electromagnetics [8]. 

Remark 2: If Vout were to be occupied by a homogeneous, linear, bianisotropic material 
different from the one occupying Vm, the process of deriving Eqs. (8) and (13) shows that 
the Post constraint would apply to that material too. 

Remark 3: The Post constraint applies at every frequency, positive, negative, or zero. 
Therefore, it also applies in the time domain. 

Remark 4: Even though the existence of T may be deduced from crystallographic symme¬ 
try [11,13] or other non-electromagnetic means, it would still be filtered out by the Maxwell 
equations. 

Remark 5: The so-called perfect electromagnetic conductor, characterized by D = T5 
and = —TE [27,28], does not have a recognizable existence as an electromagnetic medium 
because it violates the Post constraint. 


3 Surface states 

The Post constraint does not negate the various manifestations of the linear magnetoelectric 
phenomenon as enshrined in the constitutive dyadics a and P, but it does constrain the 
constitutive characterization of that phenomenon. In contrast, as is clear from Eqs. (11)-(13), 
T does not appear in the Maxwell equations in the two regions as well as in the boundary 
conditions, when the constitutive relations are written for the correct induction fields (D 
and H_) in terms of the correct primitive fields (E and 5). So, we sought an alternative 
mechanism that would preserve the Post constraint but could be invoked to interpret certain 
characteristics of experimental measurements on magnetoelectric materials. 

As in Sec. 2, let Eqs. (1) hold in Vout while the constitutive relations 

D_{r,u) = e{uj)'^{r,uj) + a{u)'^{r,uj) 1 
MSLi^) = P{^) ’ E_{r_,io) + K{io) ' B_{r,u) J 

along with Eqs. (3) hold in Vin- As boundary conditions. 


L ^ Vin I 


(15) 


let us use Eqs. (13) without 
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assuming that S is charge neutral and current neutral; i.e., 


fkiHs) X [EoutiHs,^) - = 0 

fk{r,) • u) - g{uj) • E^^ir^, u) - ^(o;) • a;)] = p,(r„ a;) 


h(r J X UoB^^t{r„ u) - (3{u) • a;) - u{u) • a;) = a;) 




(16) 

We posit that the bianisotropic material possesses surface states due to the abrupt ces¬ 
sation of its microscopic morphology, and that these surface states can be described macro- 
scopically as 


Psir^,uj) ='y{u)h{r^) ' Bi^{r^,uj) 


(17) 


where 7 is an admittance characteristic of the material. What may be observed experimen¬ 
tally with magnetoelectric materials are not the effects of a constitutive parameter that 
describes electromagnetic phenomenons in a volume, but of an admittance 7 that describes 
electromagnetic phenomenons on a surface. 

Remark 6: If Vout were to be occupied by a homogeneous, linear, bianisotropic material 
different from the one occupying Vj„, and both materials were to possess surface states, then 


Ps(r„a;) =h(rj • [7m(a;)5.„(r,, a;) - 7o„4(a;)5„„i(r,, a;)] 

= -IkilLs) X [7m(w)^i„(r„a;) - -iout{uj)E^ut{lLs^^)] 


(18) 


because the unit normal vector at G 5 pointing into Vin is —fR(Ls)- 


4 Illustrative Example 


In order to illustrate the effects of 7 , let us consider a simple boundary-value problem that 
has been put forward to indicate the measurability of [ 6 ]. Let us actually begin with a 
slightly more complicated problem [25] to be simplihed later. 

Suppose that Vout is the half space z < 0, Vin is the half space z > 0, and therefore S is 
the plane z = 0. Furthermore, let e = eJ, ^ a = 0, and /3 = 0, where the dependency 

upon u has been omitted for brevity. 

Let the primitive electromagnetic helds in Vout be 


Ml) 

Ml) 


Hy "h 


rsUy + Tp 


_ ^ 
LO 


+ 


—ttp Uy -|- a 

-rpUy + Vs 


— u ^ tq + u^k 
ko 

u ^ tq + u ^ k . 
ko 

* \ ko 

u T-o+U K, 


ex.]){iToz) 
exp(—zTo^;) \ exp(ma:) 


ko 


exp(iroz) 

exp(—iToz) [ exp(mx) 


, 


z<0. (19) 
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where and Tq = +^^0 ~ Representing a plane wave incident on 5, the 

coefficients a* and are presumed to be known. Representing the plane wave reflected 
into Vouti the coefficients and are unknown. Equations (19) satisfy the homogeneous 
counterparts of Eqs. (5). 

The primitive electromagnetic fields in Vm are given as 


E{r) 

Mir) 


ts iky d" tp 


~tp Uy T ts 


exp{iTz) exp{iKx) 
exp(irz) exp(ma;) 


z>0, ( 20 ) 


where k = A;o'\/(z^ 7 ^\/(^ 7^5 = +\/P — and the coefficients tg and tp are unknown. 

Representing the plane wave refracted into these expressions satisfy the homogeneous 
counterparts of Eqs. ( 6 ) and ( 8 ). 

The foregoing expressions were substituted in Eqs. ( 16 ) 2,4 and (17)2 to determine r^, r^, 
tg, and tp in terms of ttg and Op. The results found are as follows: 


r 


S 


ts 

tp 


{r]r - (5^)(1 + rjrSr) - {'IVofVr^r _ . 

iVr + <5r)(l + VrSr) + (7^o)^^r'^r {Vr + <5r)(l + Vr^r) + (7ho)^hr^r ’ 

{r]r + (5r)(l - r]r6r) + {-fVofVr^r _ ‘^IVoVr^r _ 

{rjr + (5r)(l + Vr^r) + {iVo^Vr^r (hr + <5r)(l + Vr^r) + {iVo^Vr^r ’ 

_ 2hr-(l +hr<^r) _ _ 2-fr]or]‘^6r _ , . 

(hr + ^r)(l + hr<5r) + (7ho)^hr<^r (hr + ^r)(l + Vr^r) + (7ho)^hr<^r ’ 

_ 2hr(hr + ^r) ___ /24i 

(hr + <5r)(l + hr<5r) + (7ho)^hr^r ^ (hr + <5,.)(1 + hr<5r) + {iVo^Vr^r ’ 


where 




(25) 


We have verified that Eqs. (21)-(24) satisfy Eqs. ( 16 ) 1,3 and (17)i; furthermore, these ex- 
presssions simplify to the standard results 


— „ Tjr—Sr 

' S - Ug , 

f — Q "^Vr 

^^Vr+Sr ’ 


— „ l—r)rSr 

' P “P l + TJrSr 

f = n _ "ilt. _ 

''P “P 1+yrSr 


(26) 


for 7 = 0 [29]. 

Cross-polarized reflection in the foregoing problem has been used to adduce the invalidity 
of the Post constraint [25], but Eqs. (21)-(24) show that both cross-polarized reflection and 
refraction are due instead to the surface states as accommodated phenomenologically through 

7- 

Let us now simplify the results to apply to the boundary-value problem put forward to 
indicate the measurability of T [ 6 ]. For that problem e = Eq and v = Uq. Then, 6r = rjr = 1 




























and Eqs. (21)-(24) simplify to: 


r. = 


aslVo-ap 

' 4 +( 7 t;o)^ ’ 


P ' '0 4+(7?7o)^ 


. _ r^2a^+apjrjo 

^ 4+( 7 ^ 0 )^ ’ 


tp 


o 2(2p—<2.s7^0 

^ 4+(7r;o)^ 


(27) 


As both half spaces z < 0 and z > 0 are vacnous, Eqs. (27) indicate that (i) reflection occurs 
solely due to non-zero 7 and (ii) the polarization of the reflected plane wave is rotated with 
respect to that of the incident plane wave. Both characteristics had been attributed to a 
non-zero 4/ [ 6 ], but that is clearly unnecessary. 


5 Concluding remarks 

We formulated the frequency-domain Maxwell equations for a homogeneous, linear, bian- 
isotropic material occupying a bounded region, as well as for the vacuous region that makes 
up the rest of space. The axionic piece T is absent in vacuum and it disappears from the 
Maxwell equations in the bianisotropic material when the induction fields ^ and ^ are re¬ 
placed by the primitive fields E and 5 by using the constitutive equations. In consequence, 
the axionic piece also vanishes from the boundary conditions prevalent at the interface of 
the two regions. The Post constraint T = 0 is thereby vindicated. 

Even if both regions were were to be occupied by dissimilar homogeneous, linear, bian¬ 
isotropic materials, the derivation process indicates that the Post constraint would apply to 
both materials. Even though the existence of T may be deduced somehow, it will be filtered 
out by the Maxwell equations. Furthermore, every biisotropic material must be Lorentz re¬ 
ciprocal and the so-called perfect electromagnetic conductor does not have a recognizable 
existence as an electromagnetic medium because it violates the Post constraint. 

We also found that characteristic effects that may be observed experimentally with slabs 
of a magnetoelectric material are not the consequences of the axionic piece of the linear 
constitutive dyadic of that material. Instead, they are the consequences of an admittance 
that describes surface states arising from the abrupt cessation of the microscopic morphology 
of the magnetoelectric material. This admittance will also impact classical analyses [30,31] 
of light scattering by topological insulators. [32] 
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